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Abstract. Our adaptive optics observations of nearby AGN at spatial resolutions as small as 
0.085" show strong evidence for recent, but no longer active, nuclear star formation. We begin 
by describing observations that highlight two contrasting methods by which gas can flow into 
the central tens of parsecs. Gas accumulation in this region will inevitably lead to a starburst, 
and we discuss the evidence for such events. We then turn to the impact of stellar evolution on 
the further inflow of gas by combining a phenomenological approach with analytical modelling 
and hydrodynamic simulations. These complementary perspectives paint a picture in which 
all the processes are ultimately regulated by the mass accretion rate into the central hundred 
parsecs, and the ensuing starburst that occurs there. The resulting supernovae delay accretion 
by generating a starburst wind, which leaves behind a clumpy interstellar medium. This provides 
O j i an ideal environment for slower stellar outflows to accrete inwards and form a dense turbulent 
disk on scales of a few parsecs. Such a scenario may resolve the discrepancy between the larger 
scale structure seen with adaptive optics and the small scale structure seen with VLTI. 
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| 1. Introduction 

Many studies of active galactic nuclei (AGN), particularly those concerned with under- 
standing the co-evolution of black holes and their host galaxies through cosmic time, rely 
on observations of quasars. The reason is simply that these are the most luminous of such 
I . objects and can be studied at all redshifts. While they form an important component of 
' the AGN population, the luminosity function (Hao et al. 2005) indicates that the num- 
ber density of lower luminosity AGN is orders of magnitude greater. This is particularly 
true in the local universe because of the effects of downsizing, which means that at low 
redshift black hole growth occurs primarily in AGN with black hole masses < 10 8 M© 
(Best et al. 2005) and Seyfert-like luminosities (Hasinger et al. 2005). As a result, with 
only a few exceptions such as Mkn231 (Davies et al. 2004), detailed studies of nearby 
" ' AGN have to focus on the lower luminosity AGN such as Seyferts. 

This has important implications. Although most quasar activity may be fuelled via 
mergers of gas rich galaxies, Seyfert nuclei typically reside in (preferentially early type) 
spiral galaxies (Ho 2008). As a result, Seyfert activity is likely to be fuelled via secular 



processes associated with disk evolution such as those outlined by|Shlosman et al. (1990) 



and Wada (2004) One obvious part of the overall scheme is the role played by bars, which 



are known to drive gas inwards from the large-scale disk. Yet despite the many attempts 
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to find a correlation between the presence of a bar and an AGN, the statistical outcome 
is, at best, for only a marginal link between these phenomena. This, perhaps, should 
not be surprising. Bars drive large amounts of gas inwards (of order 0.1 M Q yr _1 , Regan 
& Teuben 2004), working on spatial scales of kpcs and timescales of Gyrs. In contrast, 
a Seyfert nucleus requires typically < 0.01 M Q yr -1 of gas (Jogee et al. 2006), and the 
relevant temporal and spatial scales are Myr and pc. The two regimes are orders of 
magnitude different in every respect. Thus, while it is well understood how bars can 
create a gas reservoir and a starburst in the central kiloparsec (and there is plenty of 
observational evidence to support the link between bars and both of these phenomena: 
Sakomoto et al. 1999, Laurikainen et al. 2004, and Jogee et al. 2005 are a few exmaples 
among many), this is only a single possible step in the process that leads to AGN fuelling. 

It is now understood that the properties of the central kpc - the circumnuclear region 
- can play an important role in driving gas further inwards. But although observations 
have yielded many tantalising results, reaching definitive conclusions is difficult. For 



example, Martini et al. (2003) looked at the detailed circumnuclear dust morphologies of 
matched samples of active and inactive galaxies. One of their results was that AGN are 
only found in galaxies that exhibit dusty spiral patterns. But they also found that spiral 
patterns of any sort were associated equally often with active and inactive galaxies. This, 
again, is likely a result of the differing spatial and temporal scales of the phenomena. 
The aspect of time variability is one of the issues addressed by NUGA (Garcfa-Burillo et 
al. 2007), a study of the distribution and kinematics of the molecular gas in nearby low 
luminosity AGN. One of their key results was that gravitational torques can not only 
drive gas in to form a circumnuclear ring, but also tend to drive gas out from the nuclear 



region towards the ring. This led Garcfa-Burillo et al. (2005) to propose a scenario in 



which further inward flow was possible episodically if two conditions were fulfilled: the 
bar should be weakened by dynamical feedback, and the viscosity should increase due to 
a sufficiently steep density gradient in the circumnuclear ring. 

These examples show that there is still much we do not understand about the conditions 
under which gas can be transported inward from kpc scales to the nuclear region, and 
it is one of the topics addressed here: we describe two highly contrasting case studies 
where gas is being driven inwards to the central tens of parsecs. We also show evidence 
that, once there, the gas is likely to trigger a starburst. And we discuss the impact of the 
starburst on the fate of the gas, and a possible link to the sub-parsec scales. 



2. Watching Gas Flow Inwards 

In this section, we describe contrasting cases where gas can be observed as it flows in 
towards an AGN. Although we present only two examples here, analyses of non-circular 
motions and gas inflow mechanisms on similar scales have been performed for other 
galaxies as well: NGC7469 (Davies et al. 2004), NGC6951 (Storchi-Bergmann et al. 
2007), NGC4051 (Riffel et al. 2008), NGC7582 (Riffel et al. 2009), NGC4151 (Storchi- 
Bergmann et al. 2009). 

In one case, NGC 1068 (Mueller Sanchez et al. 2009), the inflow rate is prodigious, due 
to chance events, and probably cannot be sustained. In the other, NGC 1097 (Davies et 
al. 2009), the inflow is much more moderate and appears to be regulated in a way that 
is sustainable for timescales of a Gyr. Intriguingly, in neither case do we see gas flowing 
all the way down to the AGN: the inflow is terminated at spatial scales of a few to tens 
of parsecs, where it is instead expected to fuel a starburst. 
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Figure 1. Flux and velocity maps of the 1-0 S(l) H2 emission line in NGC 1097. Far and centre 
left: flux map and residual after subtracting elliptical isophotes, showing the 3 photometric arms. 
Centre and far right: velocity fiel d and residual after subtracting circular motions, showing the 
2 kinematic arms. Adapted from Davies et al. (2009) 



2.1. Gas inflow in NGC 1068 

Among the nearby AGN we have studied (Davies et al. 2007, Hicks et al. 2009), NGC 1068 
is rather unusual. The H2 emission in the central 250 pc originates in an expanding off- 
centre ring or shell with particularly bright and massive clumps around the north east 
side. Filaments of gas extend from the ring at a radius of about 30 pc to the AGN on 
both sides. Mueller Sanchez et al. (2009) has modelled the morphology and kinematics 
of the filaments, finding that the only way to simultaneously account for both constraints 
is if the filaments trace gas that is falling almost directly in towards the AGN. These 
models indicate that the infall timescale is about 1.3 Myr. One of the filaments lies across 
the front of the AGN, suggesting that inelastic collisions may allow the gas to settle on 
scales of a few parsecs. 



Estimating the mass of gas in the filaments is extremely difficult. Mueller Sanchez et al. (2009) 
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used a variety of methods to derive a mass of ~ 2 x 10 M©, although with about a factor 
of 3 uncertainty. Thus the inflow rate to the central few parsecs is of order 15 M© yr _1 . 
This is a remarkably high rate, 2-3 orders of magnitude greater than that needed to 
power the AGN itself. And because the characteristics of the gas in the central region are 
so unusual, and the inflow appears to be due to chance combination of circumstances, 
such a rapid inflow is probably unsustainable. 

2.2. Gas inflow in NGC 1097 

NGC 1097 exhibits many classical features of disk galaxies: a large scale bar, a circumnu- 
clear ring and starburst, an inner spiral, and a low luminosity AGN. Its one peculiarity 
is that the inner spiral has 3 photometric arms (Prieto et al. 2005). Integral field obser- 
vations (Davies et al. 2009) showed that this pattern is due to obscuration of stellar light 
by gas and dust in the molecular arms. Intriguingly, as seen in Fig.[T] the gas kinematics 
reveal a strong non-circular velocity residual in the form of a 2-arm spiral - as expected 
from linear theory which indicates that the projected line-of-sight velocity pattern of an 
m-arm spiral is an (m — l)-arm spiral. The properties of the arms allow one to derive 
an inflow rate along them of ~ l^Mgyr -1 . However, hydrodynamical simulations of 
nuclear spirals show that there is also significant outflow between the arms in the low 
density regions (gas from the arms overshoots the nucleus and continues as a diverg- 
ing outflow). In the simulations, the net inflow rate is a factor 20 less; applying this to 
NGC 1097 means it could be as little as 0.06 M© yr" 1 . 

Given that there is of order 10 8 M© of gas within a radius of 700 pc, such a modest 
inflow rate would require 1.8 Gyr to drain the gas reservoir in the circumnuclear region. 
If this gas is replenished by the ring, which itself is fed by the large-scale bar, then the 
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Figure 2. Relation between the characteristic age of mo st recent starburst a nd the AGN lu- 
minosity and accrection rate. This figure is adapted from Davies et al. (2007) and includes an 



additional point for NGC 2992 from IFriedrich et al. (2009) 



inflow could in principle continue even longer. Thus the inner spiral in NGC 1097 appears 
to be driving gas inwards at a sustainable rate. 

As for NGC 1068, the gas is not reaching the AGN. Instead the spiral arms appear to 
terminate in the central 10-20 pc, at which point the gas dispersion increases dramatically 
(Davies et al. 2007, Hicks et al. 2009). If the gas piles up here, it will inevitably lead to 
a starburst - and a nuclear starburst has been reported (Storchi-Bergmann et al. 2005, 
Davies et al. 2007). The estimated gas inflow rate could lead to episodic starbursts similar 
to this one every 20-150 Myr. 



3. Nuclear Starbursts 

We have seen in the previous section how, at least for these 2 cases, gas flows into the 
central few to tens of parsecs. If it is stalled here, one would naturally expect it to trigger a 



starburst. Evidence for starbursts on these scales was presented by Davies et al. (2007) 



who used near infrared integral field spectroscopy combined with adaptive optics to 
achieve high spatial resolution. Using the spectral information, these authors were able 
to separate out the non-stellar continuum (hot dust emission) associated with the AGN, 
and show that in every case the stellar continuum itself was spatially resolved with a 
typical size scale <50pc. Furthermore, they showed that this was a distinct stellar pop- 
ulation, younger and dynamically cooler than the bulge: photometrically, from an excess 
in the stellar continuum above an r 1 / 4 law fitted to 1-2" scales and extrapolated in- 
wards; and kinematically, from a drop in the stellar dispersion on the same spatial scales. 
Constraints on stellar population synthesis models provided by several independent diag- 
nostics indicated that the starburst ages were young, lying in the range 10-300 Myr. One 
important piece of evidence was the remarkably low equivalent width of the Br7 line, 
which shows unambiguously that the star formation has ceased: i.e. although still young, 
these are post-starbursts. Intriguingly, starburst models imply that during the short (per- 
haps only 10 Myr) actively star forming phase, the luminosity of the nuclear starburst 
would have been an order of magnitude greater than it is now - with a star formation 
rate per unit area comparable to that in ULIRGs. This high intensity appears to make 
sense in terms of the high gas mass surface density which, from the Kennicutt-Schmidt 
law, is expected to lead to high star formation rates. 

More pertinent to the issue of gas inflow, the authors showed that there appears to 
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be a relation between the characteristic age of the most recent episode of star formation 
and the luminosity of, or accretion rate onto, the AGN (Fig. [2]). This figure suggests that 
there is a delay of 50-100 Myr between the short duration starburst and the onset of 
accretion onto the AGN. This timescale is interesting because it corresponds to the age 
at which, for a short burst, one expects the typeii supernovae phase to be finishing. This 
topic is discussed in more detail in Section 01 which addresses the impact of stellar winds 
and supernovae on gas inflow and outflow, and how these might provide a link between 
phenomena observed with adaptive optics on ~10pc scales and those seen with infrared 
interferometric techniques on 0.1-1 pc scales. 



4. Stellar Outflows and the Torus: Linking lOpc and lpc Scales 

The data presented in |Davies et al. (2007)| and |Hicks et al. (2009) show that nuclear 
star formation is occuring on the same spatial scales, up to a few tens of parsecs, as the 
gas distribution; and that on these scales the two components have similar kinematics. 
This means the stars and gas must be mixed. It is inevitable that the star formation will 
have some impact on the gas, and this was argued from a phenomenological approach by 



Davies et al. (2007) Hydrodynamical simulations of the impact of stellar evolution, such 
as those presented in Schartmann et al. (2009) |give a clearer picture. To understand this 
in the context of Seyfert nuclei, these simulations have been scaled to match a typical 
starburst there, as described by Schartmann et al. (these proceedings): Mbh ~ 10 7 M Q , 
M s tars ~ 2 x 10 8 M Q out to 50 pc, a dispersion of a* ~ 100 km s _1 as well as rotation, and a 
mass loss rate ~ 0. 1 M Q yr^ 1 corresponding to an age of 50 Myr but which of course would 
vary with time. The simulations show that there are two, fairly well defined, regimes 
separated by a stellar outflow speed of about 250 km s -1 . Fast outflows (which could 
represent winds from OB stars, as well as type n and type I supernovae) lead to a starburst 
wind which efficiently removes the diffuse ISM. Slower outflows (which represent, for 
example, planetary nebulae or winds from AGB stars) lead to a large number of clumps 
that can accrete efficiently to smaller scales, giving rise to a dense and turbulent disk 
with a radius of 0.5-1 pc. Taken together, these imply that there is a phase in a starburst 
when accretion to smaller scales is possible and efficient. For a short starburst, this phase 
would begin after about 50 Myr, once there are no more OB stars or type n supernovae; 
and it would end after about 300 Myr once typei supernovae start to appear. 

The size scale of the central disk is similar to that of the compact 8-12 /^m structure 
in NGC 1086 (the galaxy to which the simulation is scaled) by Raban et al. (2009) us- 



ing MIDI on the VLTI; and within a factor of a few of the compact structure seen in 
Circinus by Tristram et al. (2007) using the same instrument. As already suggested by 



Davies et al. (2008) these simulations imply that stellar outflows are a link between the 



gas structure seen on scales of 10-50 pc using adaptive optics, and those seen on 0.1-1 pc 
scales using infrared interferometry. This scenario has a number of similarities to the 



simulations of Wada et al. (2009) Although the initial conditions and the distribution 
of energy input are rather different, the result is a vertically thick (up to 10 pc from the 
disk plane) distribution at large radii, and in the inner few parsec a compact turbulent 
disk. 

If one considers either of these structures to be associated with the molecular obscuring 
torus, then one must also conclude that they both are. As illustrated in Fig. [3] and 
stated by Hicks et al. (2009)] this implies that the torus comprises multiple components 
each of which fulfils different roles; that it must be a star forming torus; and, because 
the starbursts are episodic, that it must be a dynamical entity. What we see on 10- 
50 pc scales is the global structure, a diffuse and clumpy envelope where there are short 
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Figure 3. Composite figure showing how hydrodynamical simulations of stellar outflows might 
provide a link between the 10 pc and 1 pc scale structure withi n the torus. Panels are taken from 
Raban et al. (2009) | |Hicks et al. (2009) Davies et al. (2009) and simulations similar to those 



Schartmann et al. (2009)1 but scaled to a typicai Seyfert nucleus. This figure is illustrative 



and the individual panels do not necessarily refer to the same specific object. 

repeated bursts of star formation. The slow stellar winds accrete down to smaller ^Tpc 
scales where they form the more compact structure, a dense and turbulent disk. The 
evolution of this disk - in terms of further star formation and gas inflow - is addressed 
by Schartmann et al. (these proceedings). 

One final point worth noting is the similarity between these simulations and those per- 



formed for the Galactic Center, on scales about 100 times smaller, by Cuadra et al. (2006) 



and Cuadra et al. (2008) Both simulations were performed within a 0.5 pc volume in 
which the stars, their orbits, and their outflows, matched as precisely as possible the 
known configuration at that time. In the earlier version, where there were both slow 
^200 km s -1 and fast ^700 km s" 1 stellar winds, the result was that much of the ISM 
was blown out but that clumps of cooler dense gas accreted to form a compact disk that 
extended out to ~0.04pc. In the later version, the updated stellar wind speeds were all 
fast, and nearly the entire ISM was blown out. It is reassuring that qualitatively, these 



simulations yield the same result as those of Schartmann et al. (2009) even though the 



scales are drastically different, because in both cases it is the same processes (i.e. stellar 
outflows) that govern the evolution of the region. 



5. Conclusions 

We have presented a brief review of the current status of our work on the impact that 
star formation has on gas inflow to Seyfert nuclei. We are not yet able to follow gas all 
the way down to the AGN itself, and instead we have focussed on spatial scales from a 
few hundred parsecs down to the central parsec. Our main conclusions are: 

• We have seen 2 contrasting methods of gas inflow that feed a nuclear starburst. 
While showing how gas might reach the central few to tens of parsecs, this also raises 
questions about the sustainability of inflows, how far in they reach, and whether it is 
possible to define what a typical inflow mechanism or rate is. 

• There are short intense starbursts in the central tens of parsecs around Seyfert nuclei. 
In general, starbursts will be an inevitable consequence of gas inflow. 

• Stellar outflows play a key role in the evolution of the ISM on these scales. Specif- 
ically, fast outflows tend to blow out the diffuse ISM, while slow outflows can lead to 



The Impact of Nuclear Star Formation on Gas Inflow to A GN 



125 



efficient accretion of gas onto smaller scales. This may be a link between what is seen 
on 10-50 pc scales with adaptive optics, and what is observed on 0.1-1 pc scales with 
infrared interferometry. 

• The molecular obscuring torus most likely consists of multiple structures on different 
scales that fulfil different roles. Since the state of the torus is strongly affected by the 
episodic star formation that occurs within it, we need to think of it as a dynamic entity. 
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